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DEVELOPMENT OF DIRECT-INVERSE 3-D METHODS FOR APPLIED TRANSONIC 
AERODYNAMIC WING DESIGN AND ANALYSIS 


I. Introduction 

This report covers the period -from July 1, 1988 thru December 3i, 1988. The 
primary task during this period has been the comtinued development of the TAW5D 
transonic inverse wing design method with viscous interaction effects included. 

II. Personnel 

The staff associated with this project during the present reporting period 
have been Dr. Leland A. Carlson, Principal Investigator, and Robert R. Ratcliff, 
Graduate Research Assistant. The work performed during this phase of the 
project is primarily associated with the Master's Thesis of Mr. Ratcliff. It is 
anticipated that Robert will receive his degree in May 1989. 

III. Research Progress 

During the past six months, most of the research effort has been concentrated 
in the areas of program refinements, supression of divergent spanwise 
oscillations, and viscous studies. To a significant extent, these efforts have 
been summarized in an extended abstract of a paper proposed for the AIAA 7th 
Applied Aerodynamics Conference; and this extended abstract is included in this 
report as Appendix A. In addition, comments concerning these various areas are 
presented as follows: 

A. Program Refinements — 

The following improvements have been made in the TAW5D design/analysis 
program: 

(a) A simple plotting subroutine using basic Fortran commands has been 
written to allow the user to view the designed section shapes as part of the 
output file generated by TAW5D. This addition allows the user to conduct the 
design process more efficiently. 

(b) Close to optimum damping parameters and flow iterations between surface 
updates have been obtained through numerical experimentation for one design 
case including viscous interaction using Wing-A on the medium and fine grids. 
The effect on the final design of the number of iterations before the first surface 
update has also been studied. It has been found that as the design region is 
moved closer to the leading edge that more flow iterations are needed before the 
first surface update is performed in order to insure an accurate nose shape. 
Since large changes in the airfoil occur with the first surface update, a more 
accurate solution will speed up convergence. The results from this study will be 
shown in the final report. 

(c) In the calculation of the quantity tyShV)/(/»hU) on the surface of the wing, 
which is intgegrated to yield the section shapes, the idea of evaluating fihU(i,ky,k) 


using the more conveient quantity phU<i,ky-i,k), where ky is the jth grid point 
lying on the wing, has been investigated. This method assumes that the -flux ohU 
is constant in the region immediately above the wing. In order to verify this 
approach, a subroutine was written to calculate the -flux _/9hU(i,ky,k) dirrectly 
using -finite differences. Although the numbers were slightly different from the 
first approach, the effect on the final design was negligible for the cases studied 
to date. 

B. Supression of Divergent Spanwise Oscillations — 

A? discussed in Appendix A of this report, an annoying spanwise oscillation 
of the residual, and hence of the section shapes, originates at the first spanwise 
station from the root and propagates spanwise as a decaying oscillation. This 
oscillation is due in part to the compensation terms which include spanwise 
derivatives of the potential, and it becomes more significant as the wing aspect 
ratio is decreased or the sweep is increased. 

As indicated in previous progress reports, considerable effort has been 
expended over the past year in attempts to eliminate or at least mitigate this 
annoying and sometimes troublesome oscillation. During the present reporting 
period four methods have been formulated to supress this oscillation problem. 
These are as follows: 

(1) Specify the inverse boundary condition at every other spanwise station and 
linearly interpolate the inverse displacements to the stations lying in-between; 

(2) Specify the inverse boundary condition at every spanwise grid station but 
only calculate displacments at every other station and linearly interpolate the 
displacements to the stations in-between; 

(3) Smooth the slopes, (phV)/(phU), in the spanwise direction in the design region 
before the integration of the tangency condition; and 

(4) Calculate all spanwise derivatives of the potential based upon a locally 
smoothed potential function in the spanwise direction. 

Each of these methods has certain advantages in different design situations. 
For instance, methods (3) and (4) give the designer the most flexibility in that the 
desired pressure distributions can be imposed at every spanwise grid station and 
the section shapes corresponding to each grid station can be calculated relatively 
independently of the adjacent stations. On the other hand, since the first two 
methods require interpolation, the section shapes determined by these methods 
at say odd stations are directly dependent upon the shapes at even stations. 
Although with these latter methods the designer loses some flexibility, a 
smoother distribution of section thicknesses in the spanwise direction is 
obtained, which in many cases may be an asset. Of course, from a designer's 
standpoint, method (i) is the most restrictive of the four techniques; but it yields 
the smoothest designs in the spanwise direction. In addition, and perhaps more 
importantly, method (i) converges the quickest. These methods have all been 
tested on medium grids, and studies are currently in progress to test them using 
fine grids to insure that they are independent of grid size. 
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C. Viscous Studies — 


In previous studies when viscous pressures were used as input to an 
inviscidly designed wing, the resulting section shapes were "unsatisfactory" at 
the edges of the design regions in that there was a sudden increase in thickness 
corresponding to the boundary layer inherently included in the region designed 
with viscous pressures with the consequence that the resulting profile was not as 
accurately obtained as usual. In order to determine more specifically how the 
outboard and inboard regions of the wing, which were not being desgined, affected 
the design of the middle section, a study was conducted in which the entire wing 
was designed inviscidly using pressures obtained from a viscous analysis of the 
target wing, Lockheed Wing-A. The resulting wing, while being a little thicker in 
regions of large boundary layer thickness, had a smoothly varying thickness 
distribution which corresponded closely to the displacement surface of the actual 
wing. This study verified that the solution in the inboard and outboard direct 
regions of the wing had, in fact, been responsible for the shape changes at the 
edges of the design region in the case reported in a previous paper, where only 
part of the wing was designed inviscidly. 

During the course of this project, there have been many instances where even 
a viscous analysis of a wing could not be accomplished succesfully without 
diverging due to large separated regions being predicted by the boundary layer 
solver. This phenomena frequently occurred when studying Lockheed Wings B and 
C, but sometimes it even occurred with Wing A. In fact, in some cases the 
boundary layer solver could not even be executed completely due to its inability 
to find an appropriate starting solution. In some cases these problems have been 
avoided by finding an appropriate sequence of flow field iterations and boundary 
layer updates or by going to a finer grid from the one currently being utilized. 
However, in many cases a few degrees of dihedral had to be added to the wing to 
allevite this premature separation problem. At the present time, it appears that 
this boundary layer solver problem needs further investigation. 

IV. Current and Future Efforts 

Obviously, the primary effort during the next reporting period will be the 
completion of Robert Ratcliff's thesis and the preparation of the project final 
report. However, depending upon the time and financial resources available, the 
following items may be investigated: 

(1) It has been discovered that the skewness of the grid lines leaving the wing 
tip in the spanwise direction has a large effect on the accuracy of the sections 
designed in the tip region. As a consequence, some check for grid skewness or a 
method to minimize the effect of skewness on the final solution is needed. If time 
permits, appropriate sections will be added to the program. 

(2) The logic necessary to design a segment of the wing which starts aft of the 
leading edge and terminates prior to the trailing edge is currently being studied 
and hopefully will be included in the final version of the code. If at all possible, 
test cases demonstrating this option will be shown in the final report. 
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(3) Since the spanwise stability problem is now better understood, the trailing 
edge of the designed air-foil sections may be allowed to -float i-f the correct 
method to supress the spanwise oscillation is utilized. This approach would 
allow an additional degree of -freedom in the design solution procedure and may 
permit the solution o-f "difficult" design problems. This conjecture will be tested 
by including logic which allows the trailing edge to "float" and performing 
appropriate test cases using the previously discussed spanwise oscillation 
suppression methods. In these tests, the accuracy of the designed sectional 
thickness profiles and the sectional angle of attacks will be ascertained. 

(4) As mentioned in previous reports, two different methods are used in FLO30 to 
define the potentials at ghost points in the fuselage adjacent to the root section 
of the wing. One approach uses the flow tangency condition, while the other 
linearly extrapolates the potentials in the spanwise direction to the ghost point 
loction. Since these potentials are used in the inverse boundary conditions, their 
method of determination may have an effect on the design of the wing, especially 
at the root. Hopefully, resources will be available to permit the investigation of 
the influence of these ghost potentials. 

(5) Currently, the only time the spanwise derivatives of the potential are based 
upon a spanwise smoothed potential function, is in the calculation of the residual 
immediately before the integration of the flow tangency condition or surface 
update. Note that this procedure is only used when the fourth method of 
spanwise oscillation suppression is utilized. It is planned to investigate the 
effect on the final design of carrying out this smoothing approach every time the 
residual is calculated. 


V. Grant Monitor 

The NASA Technical Monitor for this project is Richard L. Campbell, Applied 
Aerodynamics Group, Transonic Aerodynamics Branch, Applied Aerodynamics 
Division, NA.SA. Langley, (I think — I have not kept up with recent organizational 
changes). 
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APPENDIX A 



A DIRECT- INVERSE TRANSONIC WING DESIGN METHOD 


IN CURVILINEAR COORDINATES INCLUDING VISCOUS-INTERACTION 


An Extended Abstract by: 

Robert R. Ratcliff and Leland A. Carlson 
Aerospace Engineering Department 
Texas A&M University 
College Station, Texas 77843 
Phone : (409) 845 - 7541 


December 9, 1988 


Summary 

Progress in the direct-inverse wing design method in curvilinear coordinates 
has been made. This includes the remedying of a spanwise oscillation problem and 
the assessment of grid skewness, viscous interaction, and the initial airfoil section on 
the final design. It was found that : in response to the spanwise oscillation problem 
that designing at every other spanwise station produced the best results for the cases 
presented, a smoothly varjdng grid is especially needed for the accurate design at the 
wing tip, boundary layer displacement thicknesses must be included in a successful 
wing design, the design of high and medium aspect-ratio wings is possible with this 
code, and that the final airfoil section designed is fairly independent of the initial 
section. 


1 
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A DIRECT-INVERSE TRANSONIC WING-DESIGN METHOD 
IN CURVILINEAR COORDINATES INCLUDING VISCOUS-INTERACTION 

Introduction 

With the advent of efficient numerical schemes that accurately model the ir- 
rotational transonic flow about complex configurations such as wing-bodies and 
computers with ever increasing memory capacities and computational speeds nec- 
essary to execute these schemes in a reasonable amount of time, the efficient design 
of wings for transonic flight is quickly becoming a reality. Although transonic poten- 
tial schemes combined with integral boundary layer solvers may not perfectly model 
the real flow field as closely as Euler or Navier Stokes Schemes, they can signifi- 
cantly reduce the exhorbitant wind tunnel costs and time expenditures associated 
with transonic wing design. 

Many methods ranging from optimization techniques 1-2 to various inverse 
methods have been formulated using potential solvers to design wings in transonic 
flight 3-9 . One such method, which has been under development at Texas A&M 
University for the last several years, is the direct-inverse transonic wing design 
method. In this method the airfoil sections making up the wing are created by 
specifying desired pressure distributions over all or part of the w r ing aft of the leading 
edge, solving via finite-difference or finite- volume techniques the mixed Neumann 
and Dirichlet boundary value problem associated with the full potential equation 
for compressible flow, and then integrating the flow boundary condition at each 
spanwise station in the design region. The pressure distributions can be selected by 

Journal model is AIAA Journal of Aircraft. 
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the experienced engineer to have such desirable characteristics as mild or nonexistent 
shocks, a slowly increasing adverse pressure gradient, a center of pressure giving a 
desirable pitching moment or an efficient spanwise loading. The designer may also 
use wind tunnel tests of successful airfoils as an aid in picking a desirable pressure 
distribution. 

The direct-inverse technique has been successfully used in a stretched and 
sheared cartesian system 6-7 and more recently in a curvilinear system. This paper 
presents progress in the latter. It will include a brief description of the analysis 
and design methods, techniques used to suppress a spanwise oscillation problem 
resulting from the interaction of the design method with the potential solver, and 
it will present a series of test cases that will reveal the lack of dependency of the 
design on the initial airfoil section, the importance of including viscous effects in 
wing design, and any constraints due to aspect ratio and sweep of the wings and 
grid skewness in the spanwise direction. 

Background 

The base Fortran program which was modified intially by Gaily 9 for inverse 
wing design, is TAWFIVE 10 . This program not only has the capability of computing 
the potential field about a fairly general wing and fuselage combination but also 
includes a three dimensional integral boundary layer scheme to provide the necessary 
viscous effects, including the boundary layer displacement thickness, wake curvature 
and thickness. 

The inviscid numerical scheme is based upon Jameson’s and Caughey’s con- 
servative, finite- volume, full potential flow solver, FLO30. Computations are per- 
formed on a body-fitted, sheared, parabolic, wind-tunnel type coordinate system. 
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The theory behind this will be briefly discussed, but the reader is encouraged to 
refer to references (11-14) for a more in depth study. 

FLO30 solves the compressible potential equation in conservative form written 
in curvilinear coordinates : 


( P hU) ( + ( phV) v + (phW) c = 0 


( 1 ) 


where the physical velocities, u,v,w are related to the gradient of the reduced 
potential function, (j > , by : 

-1 / (Cos {a)' 


u 

V | = 

w 


H J 


<t>v — t- Sin (a) 

Ml V 0 


( 2 ) 


and the contravariant velocities U, F, W are related to the physical velocities by : 

( 3 ) 


U \ u 

V = H -1 v 
W ) \w 


where H is the transformation matrix defined as: 


H 


x i 

x v 

X C 

Vi 

Vv 

VC 

H 

z v 

Z C 


with 


h =|H| 


( 4 ) 


The local air density normalized by the freestream density can be derived in terms 
of the local speed of sound using the energy equation and isentropic relations to 
yield : 

/ „ \ 5 

( 5 ) 


P = 


a 

Q-cc 


where the local speed of sound normalized by the freestream speed, q^, can be 
calculated in terms of the speed of sound at stagnation using : 

2 


2 2 

a = a 0 - 


Qoo 


7-1 


( 6 ) 
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In a finite-volume scheme, the discrete nature of the finite difference model is 
considered from the onset and its formulation is directly analagous to the method 
in which the original compressible continuity equation was derived. The FLO30 
finite-volume scheme does this using a staggered box approach. This method used 
to discretize the differential equation is best explained with the aid of the 2-D cells 
in figure 1 and the 2-D, incompressible version of equation 1 written in cartesian 
coordinates : 

(U) X +(V) y =0 (7) 


The staggered boxes make up sets of interlocking primarj r and secondary cells. The 
values of the potential function at grid points, which are located at the corners of 
the four primary cells shown, are used to calculate the fluxes, U and V, at the 
centers of each primary cell (which are also the corners of the secondary cell) in the 
following manner : 


U = <f>x = V-y&x4> 


V = <f>Y = 


( 8 ) 


where fi and 6 are averaging and differentiating operators respectively and are 
defined as : 

f* x f = 2 + ^ 

— /i+i,; — fi — \ ,j 

where it is assumed that AX = 1 These fluxes , U and V, calculated at each of 
the primary cell centers are then averaged in the Y direction and the X direction, 
respectively, and the net flux into the cube is obtained resulting in: 


fiY&x (U) + nx&Y (JO = 0 


which is the discretized version of equation 7. 


( 10 ) 


II 
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Jameson found that lumping the fluxes at the primary cell centers led to an 
uncoupling of the solution between adjacent grid points. Therefore, he added some 
compensation terms which basically extrapolate the fluxes from the corners of the 
secondary cell to a distance, e, towards the midpoint of each secondary cell face. 
For example, given an e = .25 the flux, U, at the corresponding grid location 
(*'+ U+ J) is : 


U- 


i 1 = U-,1 

1 i 2 ' 4 * ' 2 1 


i+\ 


- .25 | ^ 

av 


*+jj+ 2 


( 11 ) 


where : 




*+5>i+5 


( 12 ) 


When all the fluxes are extrapolated in this manner and included in equation 10, it 
leads to : 


PYY&Xx4> + PXX^Yy4> — -255 X XYY <f> = 0 


(13) 


If the same procedure is extended to a three dimensional fiowfield in curvilinear 
coordinates, the discretized full potential equation becomes : 

( phU + P) + (phV + Q) + (phW + R) 

- e (^H&znQi-n + C + HhtQtt = ® 

where the Q’ s are the compensation terms defined by Jameson as : 

Qtv = (4e + ^v)^c 
QvC ~ (Af) + A{) 

(15) 

Q cf = (^c + p- v b^<t> 

Qa = + A v + A ( ) s^4> 

And, A^, A v , A( are influence coefficients which compensate for the dependence of p 
on 4>ii 4>7 it and Pi Qi an< l P are Jameson’s upwinding terms which desymmetrize 
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the scheme in supersonic zones and exclude unreal, discontinous expansions from 
the solution by providing an artificial viscosity. This equation, solved via SLOR, is 
of course a direct statement of the conversation of mass and should tend to zero as 
the solution converges unless mass is being created within the secondary cell. 

As mentioned earlier, a body-fitted, wind tunnel-type grid is used in FL030 
which is shown in figure 2. The grid shown is the coarsest mesh and has 40 x 
6x8 grid points in the X ', F, and Z directions, respectively. With this grid, 
the wing becomes a constant Y surface, and each cylindrical looking shell is a 
constant Z surface. Constant X lines can be seen running spanwise on the wing at 
a constant chord fraction from the leading edge. Notice also, due to the conformal 
transformation used, that constant X lines are packed close to the leading edge of 
the wing. This becomes to be an attractive feature when designing airfoil sections 
using the direct-inverse approach. Moreover, constant Z lines are spaced evenly on 
the wing and will, on the finest mesh, give the designer up to 21 spanwise stations 
where at the pressure distributions can be specified. 

The tangency boundary conditions are simply enforced by reflecting the fluxes 
above the surface to ghost points below the wing, fuselage or symmetry plane such 
that the out of plane component of the flux is zero at the surface. 


Inverse Design Method 

As stated earlier, in the direct-inverse method a pressure boundary condition 
is enforced rather than flow tangency aft of the leading edge of the sections of the 
wing which are to be redesigned. After Gaily 9 ’ 15 , the input pressure coefficient can 
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be written in terms of the Mach number, M , and the dynamic pressure, q, as: 


Cp 7 M\ 


1 + lzl M '(i- 9 


9 2 oo 


7-1 


- 1 


(16) 


where q 2 = (u 2 -f v 2 + w 2 ) q 2 


Solving for u in equation 2 and equation 16 and then equating the two results 


in : 


Jn <f>( + = 


1 - 


(7-1 )M 2 c 


( 1 + 7Ma ” Cf ] " 


7-1 

' - 1 


1 + («) + (5)' 


— COSO (17) 


, -1 


where Jij are the elements of ( H T ) 

A potential, 4>i,j,ki can he found in terms of the pressure coefficient by 
expanding about the grid point location (t+|,j,A:), and then using central 
differences in the f and £ direction and second order backward differences in the 
normal direction, 77 , yielding : 

^ i,ky,k ~ "T : 3 r 1 *^ 11 ^ i-l,ky,k 

i'll + 4 "12 l 

— J 12 — 4 {4^ i,ky-\,k + i—l ,ky—l ,k) 

+ < f >n i,ky-2,k + ^ n t-l,fcy-2,fc / 4 

*^13 i,fcy,fc 4 -l ” 1 ” ^ i — l,fcj/,fc+l ^ i,ky,k — i $ i — l,ky,k—l) 


where i * 1 is the right hand side of equation 17 and j — ky on the 

wing surface. 

Since the grid is boundary conforming, the wing sections in the design region 
must be updated every so often by integrating the flow tangency condition written 


14 
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in curvilinear coordinates: 


dp\ __V W f dr) \ 

A, V ~Tj \T() iitlhk 


( 19 ) 


Assuming that the spanwise terms are relatively small, this equation reduces to: 


dr i 


i,ky,k 


i,ky,k 


( 20 ) 


(This assumption was verified by solving equation 19 iteratively using constant 
values for the ratios of the velocities obtained at the present time level; the spanwise 
terms were at least two orders of magnitude smaller than the chordwise terms prior 
to the creation of the new grid, when of course all the derivatives of 77 with respect 
to the spanwise direction, (, or the chordwise direction, f, on the wing are zero.) 

V can be obtained most accurately from the residual expression. Using the 
predefined averaging and differencing operators : 


P’1 (P hV kk,,t = I ( ( P hV ),,k,+ i,k + 


*y (phV) iM , k = ((phV) iM _ iit - {phV) iM+iti 

and combining, yields : 


6y iphV) ilvl = 2(phV) i tl _, k - 2,iy (phV) iMik 


( 21 ) 

( 22 ) 


(23) 


Substituting this into the residual expression, equation 14, and solving for the out 
of plane flux, phV, on the wing surface yields : 


2 Ptv< (phV)i,ky,k = P V (S X ( phU) iM k + 2p i<: (phV ijky _ i )fc ) + 

{phWi t k y ,k) + compensation and upwinding terms 


( 24 ) 
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Now, since at convergence the flow should also be tangent to the designed 
surface, the tangency condition is enforced in the residual expression, equation 14, 
by setting : 


(P hV ) i, = ~(P hV )iM-\,k 


(25) 


Since the resulting expression is identical to equation 24, the expression for the 
normal flux becomes : 

Residual 

P'tvC = T7j 



Note that since the residual is not zero in the design region due to the inverse 
boundary condition, this expression reveals that there will be an ejection of fluid 
from the boundary. No attempt was made to account for this addition of mass 
into the flowfield, since at convergence it would be negligible. Upon substitution of 
equation 26 into equation 20 and using the cell averaged flux phU on the surface 
this becomes : 


dr] Residual 

8( 2 P(<i( (pMJ) 


(27) 


The displacements normal to the surface are obtained by integrating from the 
beginning of the inverse region to the trailing edge at each spanwise station using 
the trapezoidal rule. (Higher order integration schemes were tried but had little 
effect on the final answer, except on a fairly coarse grid in regions of high curvature 
like the cove region of a supercritical airfoil.) Assuming that the grid line leaving 
the wing is normal to the wing, this displacement, Ar], is then converted from 
computational to physical units by scaling it by the transformation metrics such 
that : 


Al = Ar)\ 


lb 


dx 7 dy 7 

M 


(28) 
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After subtracting the boundary layer displacement thickness from the inverse 
displacements, A Z, which have been linearly interpolated to the user defined input 
stations, the resulting displacements are added to the initial airfoil sections yielding 
the new wing surface for the current time level. 

Many times, the trailing edge thickness may be too large if the leading edge 
curvature is too small or may be ‘fish-tailed’ if the leading edge curvature is too large. 
These undesirable situations are remedied by a procedure called relofting where the 
designed surface is rotated about the leading edge to meet a specified trailing edge 
ordinate. To illustrate the previous procedures, a history of the convergence of a 
typical design is shown at every seven surface updates in figure 3. 

Spanwise Oscillations 

An annoying divergent spanwise oscillation problem has surfaced when design- 
ing a wing which required extensive relofting, especially when the initial section was 
thinner than the target. This led to sections which were too thick or too thin at 
adjacent constant Z grid stations (see figure 4). This problem was more pronounced 
when the sweep was increased or the aspect ratio was decreased. After many failed 
attempts at remedying this by reformulating the inverse boundary condition, atten- 
tion was directed towards the residual and the terms composing it. The residual was 
broken into its components and plotted after each surface update. Sample plots for 
a divergent case are shown at three different time levels in figure 5. As can be seen, 
the compensation terms which include spanwise derivatives of <f> at first are very 
small compared to the rest of the terms, but later tend to dominate and amplify the 
oscillation. This oscillation starts at the direct-inverse interface or, in other words, 
at the first spanwise station from the root in the design region and propagates 
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spanwise as a damped oscillation with a period of two grid spacings. Presently, it 
is thought that the mismatch in the potentials at the direct-inverse interface in the 
spanwise direction is being amplified by those compensation terms which include 
spanwise derivatives of the potential function. The residual is then undershot and 
overshot on alternating spanwise stations. This oscillation is further magnified by 
relofting, which creates a section that is too thin when the slopes defined in equa- 
tion 27, which of course are directly proportional to the Residual, are too large and 
vice-versa. Since more or less fluid has to be ejected from the section thats too 
thin or thick respectively to give the streamline approximately corresponding to 
the correct target section, the potential field at each design station is taken further 
away from the adjacent fields by the inverse boundary condition which forces an 
even further undershoot or overshoot of the residual. It should be noted that this 
problem is not soley due to the implementation of the direct-inverse technique since 
this oscillation has not been observed with the ZEBRA design code, but seems to 
be unique to the coupling of the method with the analysis code, FLO30. 

After exploring many alternatives to counter this problem, four methods have 
been discovered to damp out the spanwise oscillation. 

A) Specify the inverse boundary condition at every other spanwise station and 
linearly interpolate the inverse displacements to the stations lying in between. 

B) Specify the inverse boundary condition at every station, but again only 
calculate displacements at every other station and linearly interpolate the 
displacements to the stations in between. 

C) Smooth the slopes, in the spanwise direction in the design region. (Note: 
smoothing the integrated slopes (i.e the displacements) did not suppress the 
oscillation, but only slowed the rate of divergence.) 
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D) 


Calculate all spanwise derivatives of the potential based upon a smoothed 
potential function in the spanwise direction as follows : 


(*c) 


»>?>*+ 5 


( + 4 > i,j,k + 1 4 > i,j,k-l) 

4.0 


(29) 


Preliminary results are shown in figures 6-9 for the four different approaches. 
All four were designed inviscidly on a medium grid (60x12x16) using Lockheed 
Wing-A at an angle of attack of two degrees and a Mach number of .8. A NACA 
0012 section was used as the initial section in the entire design region which stretched 
from 30-70% semispan and began 5% aft of the leading edge and extended to the 
trailing edge. Using the RMS of the error between the target sectiona and the 
section designed as a measure of goodness, the technique A and D produced the 
best results for this case in the interior as well as at the edges of the design region, 
while for the same number of flowfield iterations, the technique C produced the 
most unsatisfactory results when compared to the target sections. The effect of each 
approach on the Residual can be seen by refering to figures 10-13. The smoothing 
approaches work well when designing in the interior of the wing, but they do not 
give satisfactory results at the root and tip of the wing where smoothing the quickly 
varying potential functions leads to large errors in the residual. 


Spanwise Grid Skewness 

Recently it was discovered, quite inadvertently, that the skewness of the 
constant X grid lines leaving the tip of the wing (See figure 14 ) had a dramatic 
effect on the design of the sections near the wing tip. If the grid was significantly 
skewed it was impossible to obtain the correct airfoil shapes in the tip region as can 
be seen in figure 15. This is thought to be due to the large errors associated with 
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a skewed grid which are revealed in the pressure distributions (figure 16). These 
errors in the solution seem to have caused the shock location to move further aft. 
Although the skeweness of the grid was quite extreme in this case, these results 
affirm the need for smoothly varying grids in wing design, at least in the spanwise 
direction. 

Boundary Layer and Wake Effects 

It was of interest to discover how significant the various viscous effects were 
in the design of transonic wings 16 . An input pressure distribution was obtained by 
analyzing Lockheed wing-A using full viscous effects ; these included boundary layer 
displacement thickness, wake thickness and wake curvature. This was considered to 
be the type of real pressure distribution which would be available to the designer. 
This pressure distribution was then used in three cases. In the first case, the wing 
was designed inviscidly. In the second case, the wing was designed without the wake 
options but included the boundary displacement thickness effects. And in the third 
case, all viscous effect were used in the design of the wing. The design region for 
all three cases extended from 30-70% semispan and began 10% aft of the leading 
edge of the airfoil, but the inverse boundary condition was only enforced at the 30, 
50 and 70% semispan station and the displacements were linearly interpolated to 
the stations delimitted. The initial airfoil section at 50% semispan was formed by 
thinning the supercritical target section by 6% and removing the cove region. The 
initial sections at the edges of the design region were the same as the target sections, 
while the remaining sections delimitted were obtained through linear interpolation. 

The results for these cases are presented in figures 17. Neglecting wake effects 
seems to have a small effect on the resulting design. The sections are little thicker 
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than the sections designed with full viscous effects. This is probably because wake 
curvature effects reduce the lift of the wing resulting in smaller boundary layer 
displacement thicknesses than those obtained without wake effects. Because of this, 
the section designed without wake effects was relofted less, leading to a thicker 
section. On the other hand, the wing sections designed inviscidly are profoundly 
different at 30 and 70% semispan, but only slightly different at 50% semispan. The 
large differences at the inboard and outboard stations are due to the influence of the 
inviscid pressures outside the design region which probably cause a large mismatch 
in the potentials at the design-analysis interface. And the remarkable agreement in 
the middle of the design region, except in the cove region where the boundary layer 
is thick, is due to the influence of the viscous boundary condition at the edges of 
the design region. (This was verified by designing the entire wing inviscidly using 
viscous pressure distributions. This case led to airfoil sections which varied smoothly 
in the spanwise direction. These results will be shown in the final paper) After the 
wings were designed, all three were then analyzed with full viscous effects to assess 
the significance of the changes made to the wing on the pressure distributions and 
to see how well these pressures matched the target pressures. Knowing that the 
wing designed with full viscous effects is correct, it is quite obvious from figure 18 
and Table 1 that the wing designed inviscidly is quite unsatisfactory. The shock is 
not far enough aft and the lift produced is sometimes 20% smaller than that desired. 

Initial profile effects 

One of the disadvantages of the inverse method is that a-priori knowledge about 
the correct shape of the leading edge must be known to achieve the correct airfoil 
shape. It was thought that because FLO30’s grid package packs grid lines close 
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to the leading edge of the airfoil, that the design could be started quite close to 
the leading edge thus relieving the designer of the burden of choosing the correct 
nose shape. Two test cases were conducted to investigate the dependence of the 
final design on the initial airfoil section. Both used Lockheed wing-A at the same 
conditions mentioned earlier for the viscous study. For the first case, the initial 
airfoils were the same as those in the viscous study. These airfoils all had leading 
edges which were in the same family as the target section. The design was started 
10% aft of the leading edge. In the second case, NACA 0012 sections were used at 
all design stations; the leading edge of these sections were not in the same family as 
the target airfoil sections. For this case, the pressure boundary condition began 5% 
aft of the leading edge. Referring to figures 19 and 20, it can be said that although 
slightly better results were obtained for the first case, the airfoils designed were 
quite insensitive to the initial section. 

Aspect Ratio and Sweep Effects 

Three cases will be shown in order to roughly delimit the range of aspect ratios 
and leading edge sweep angles. These include Lockheed Wing A and B and will 
hopefully include Wing C in the final paper. These wings have aspect ratios of 8, 
3.8, and 2.6, leading edge sweep angles of 27,35, and 45 degrees and taper ratios 
of .4, .4, and .3 respectively. The target pressure distributions were obtained by 
a direct analysis of the target wings in a fully viscous environment. The initial 
section for Wing-A was a NACA 0012, while a NACA 0006 was used for Wing-B. 
The design region for Wing-B extended from 10-100% semispan and began 5% aft 
of the leading edge. Both wings were designed using all available viscous options. 
Fine grid results for Wing-B are shown in figure 21, while medium grid results for 
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Wing-A are shown in figure 22. As can be seen the designed and target sections 
for both Wings are in excellent agreement in the interior of the design region and 
closely match at the edges of the design region. 

Conclusions 

Progress in the direct-inverse wing design method in curvilinear coordinates 
has been made. This included the remedying of a spanwise oscillation problem and 
the assessment of grid skewness, viscous interaction, and the initial airfoil section on 
the final design. It was found that : in response to the spanwise oscillation problem 
that designing at every other spanwise station produced the best results for the cases 
presented, a smoothly varying grid is especially needed for the accurate design at the 
wing tip, boundary layer displacement thicknesses must be included in a successful 
wing design, the design of high and medium aspect ratio wings is possible with this 
code, and that the final airfoil section designed is fairly independent of the initial 
section. 
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Figure 7a-7d. Airfoil Sections Designed Using Method B 
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figure 14a-14b. Comparison of a Fairly Unskewed 
and Skewed Grid at the Wing Surface 
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Figure 16a-16b. Comparison of Pressures Obtained From 
an Analysis of Wing- A Using the Grids 
in Figure 14 
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Figure 17b 
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Figure 18a-18e. Comparison of Pressures Obtained 
by a Viscous Analysis of the Wings Designed Using 
Different Interaction Assumptions 
Lockheed Wing-A, Mach - .8, AOA = 2°, Re = 2-lxl 
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Table 1. Comparison of the Wing Lift Coefficient 
for the Wings designed using different 
Viscous Interaction Assumptions 
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Figure 19a-] 9c. Comparison of the Designed Section 
with the Target Section and the First Type of Initial Section 
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Figure 21a-21c. Comparison of the Designed Section 
with the Target Section and the NACA 0012 Initial Se- tion 
I'sintr a High Aspect Ra‘io in<r 
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Figure 22a-22c. Comparison of the Designed Section 
with the Target Section and the NACA 0006 Initial Section 
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